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Abstract 

We consider a parameter estimation problem to determine the viscos- 
ity v of a stochastically perturbed 2D Navier-Stokes system. We derive 
several different classes of estimators based on the first N Fourier modes 
of a single sample path observed on a finite time interval. We study the 
consistency and asymptotic normality of these estimators. Our analysis 
treats strong, pathwise solutions for both the periodic and bounded do- 
main cases in the presence of an additive white (in time) noise. 

Keywords: Parameter Estimation, Inverse Problems, Nonlinear Stochastic Partial 
Differential Equations, Navier-Stokes Equations, Maximum Likelihood Estimators, 
Stochastic Evolution Equations, Estimation of Viscosity. 

MSC2010: 60H15, 35Q30, 65L09. 

1 Introduction 

The theory of stochastic partial differential equations (SPDEs) is a rapidly de- 
veloping field of pure and applied mathematics. These equations are used to 
describe the evolution of dynamical systems in the presence of persistent spatial- 
temporal uncertainties. When considering nonlinear processes one encounters 
many new, fundamental and mathematically challenging problems for SPDEs, 
with important applications in physics and applied sciences. 

While the general form of a particular SPDE is commonly derived from 
the fundamental properties of the underlying processes under study, frequently 
parameters arise in the formulation which need to be specified or determined 
on the basis of some sort of empirical observation. In such situations, the so 
called problem of "parameter estimation arises naturally: under the assumption 
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that a phenomenon of interest follows the dynamics of an SPDE, and given 
that some realizations of this process are measured, we wish to find the unknown 
parameters appearing in the model, such that the equations fit or predict as much 
as possible the observed data. 

Actually, the development of methods to estimate parameters appearing in 
a model serve practical considerations for two reasons. On the one hand we may 
be confident in the model, but have an incomplete knowledge of the physical 
parameters appearing therein. An "estimator" of the true parameter therefore 
provides a means to measure these unknowns. On the other hand, we may 
already possess accurate knowledge of the physical quantities involved in the 
model, but lack confidence in the validity of the underlying model. In this 
situation finding an "estimator" will be the first step in testing and validating 
the model. 

Since the solution of an SPDE is a random variable, this inverse problem of 
finding the true parameters is treated by methods from stochastic analysis and 
statistics. In this work we will follow a continuous time approach and assume 
that the solution U = U u (t,u) of the SPDE is observed for every time t over 
an interval [0, T]. We note that different types of methods and approaches are 
used to study inverse problems for deterministic PDEs, and we refer the reader 
to [20, 21] and references therein. 

A core notion in the theory of statistical inference for stochastic processes 
is the so called 'regularity' of the family of probability measures associated to 
the set of possible values O of the parameter of interest v. Note that v could 
be a vector in general. Let H be the function space where the solution evolves 
and for each v £ O denote by P^, the probability measures on C([0,T];H) 
generated by the solution U v (t), < t < T. We say that a model is 'regular' 
if any two probability measures from the family {P^, v G O}, arc mutually 
absolutely continuous. On the other hand the model is said to be 'singular' if 
these measures are mutually singular. 

For regular models one approach to the parameter estimation problem is 
to consider the Maximum Likelihood Estimator (MLE) v of v. This type of 
estimator is obtained by fixing a reference value vq and then maximizing the 
Radon- Nikodym derivative or Likelihood Ratio dP^/dP^ with respect to v. 
Usually v ^ v and the problem is to study the convergence of these estimators 
to the true parameter as more information arrives (for example as time passes 
or by decreasing the amplitude of the noise). In contrast, each singular model 
requires an individual approach, and usually the true parameter can be found 
exactly, without any limiting procedure (at least if the solution is observed 
continuously.) 

Statistical inference for finite dimensional systems of stochastic differential 
equations (SDEs) have been studied widely and provide instructive examples of 
both 'regular' and 'singular' problems. Typically estimating the drift coefficient 
for an SDE is a regular problem which may be treated with an MLE. Here 
the likelihood ratio can be determined by Girsanov type theorems. By contrast, 
estimating the diffusion coefficient is a singular problem and in this case one can 
find the diffusion coefficient by measuring the quadratic variation of the process. 
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In general there exist necessary and sufficient conditions for the regularity for 
(finite dimensional) SDEs. See the monographs [22], [25], and references therein 
for a comprehensive treatment. 

It turns out that the parameter estimation problem for infinite dimensional 
systems (SPDEs) is, in many cases, a singular problem where one can find the 
parameter v "exactly" on any finite interval of time. In particular this has been 
shown in the case of linear stochastic parabolic equations with the parameter 
of interest in the drift appearing next to the highest order differential oper- 
ator. Note that this is in direct contrast to most of the corresponding finite 
dimensional processes where one has to observe a sample path over an infinite 
time horizon or to decrease the amplitude of the noise term in order to get 
similar results. One of the first significant works in the theory of statistical 
inference for SPDEs that explorers this singularity is [17]. The idea in this work 
is to approximate the original singular problem by a sequence of regular prob- 
lems for which MLEs exist. This approximation is carried out by considering 
Galerkin-typc projections of the solution onto a finite-dimensional space where 
the estimation problem becomes regular. They prove that as the dimension of 
the projection increases the corresponding MLEs converge to the true param- 
eter. In [18, 19, 27, 28], the problem has been extended to a general class of 
linear parabolic SPDE driven by additive noise and the convergence of the esti- 
mators has been classified in terms of the order of the corresponding differential 
operators. For recent developments and other types of inference problems for 
linear SPDEs see the survey paper [26] and containing references. 

While the linear theory has been extensively studied in the framework de- 
scribed above it seems that, to the best of our knowledge, no similar results 
have been established for nonlinear SPDEs. We therefore embark in this and 
concurrent work [5] on a study of parameter estimation problems for certain 
fundamental nonlinear SPDEs from fluid dynamics. 

Note that for the linear case, key properties such as efficiency and asymptotic 
normality of the estimators, are proven by making essential use of the exact long 
time behavior of the moments of the Fourier coefficients of the solutions. In the 
case of nonlinear equations, for example stochastic equations from mathematical 
fluid dynamics, the problem is much more delicate, due to the (highly nontrivial) 
coupling of the Fourier modes. 

From the point of view of applications this work is motivated in particular by 
recent developments in the area of Geophysical Fluid Dynamics (GFD) where 
the theory of SPDEs is now playing an important role. See, for example, [31, 
32, 33, 11, 16, 15, 9]. For this developing field, novel 'inverse' methods are 
clearly needed. While the problems we consider initially are toy models in 
comparison to large scale circulation models such as the Primitive Equations, 
we are optimistic that the methods and insights developed for simple nonlinear 
SPDEs will eventually serve the wider goal of extending our understanding to 
a more physically realistic setting. 

In this work we consider the 2D Navier-Stokes equations forced with an 
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additive white noise: 

dU + ((U ■ V)U - vAU + VP)dt = adW, (1.1a) 

V-U = 0, (1.1b) 

17(0) - U , (1.1c) 



which describe the flow of a viscous, incompressible fluid. Here U = (Ui, U2) and 
P respectively represent the velocity field and the pressure. The coefficient v > 
corresponds to the kinematic viscosity of the fluid, and it will be the parameter 
of interest. The goal of our analysis will be to find a suitable estimator v = v{uj) 
which is a functional of a single sample path U(ui) observed over a finite and 
fixed time interval [0, T]. 

We assume that the governing equations (1.1) evolve over a domain T>. 
Throughout this work we will consider two possible boundary conditions. On 
the one hand we may suppose that the flow occurs over all of K 2 , take T> = 
[— L/2, L/2] 2 for some L > and prescribe the periodic boundary condition: 

U(x + Le j ,t) = U(x,t), for all x £ R 2 ,t > 0; / t/(x)dx = 0. 1 (1.2) 

We also consider the case when I? is a bounded subset of R 2 with a smooth 
boundary dT> and assume the Dirichlet (no slip) boundary condition: 

U(x,t) =0forallxe<9£>,t>0. (1.3) 

The stochastic forcing we consider is an additive space-time noise colored in 
space. Formally, we may write 

adW = J2 x r^kdW k , (1.4) 

k 

where are the eigcnfunctions of the Stokes operator, represent the asso- 
ciated eigenvalues, and Wk, k > 1, are one dimensional independent Brownian 
motions. We assume that 7 is a real parameter greater than 1 which guarantees 
some spatial smoothness in the forcing. We may also formally derive (see e.g. 
[8]) the space-time correlation structure of the noise term 

E(adW{yi, t)adW(y, a)) = Kfa y)8 t . s , 

where /v(x,y) = £ fc >! A^ 2 ^ fe (x)$ fc (y). 

We should mention that the Stochastic Navier-Stokes equations in both 2 
and 3 dimensions and under much more general stochastic forcing conditions 
have been extensively studied. See, for instance, [2, 3, 4, 7, 12, 13, 30] and 
containing references. 

Since the parameter of interest v appears next to the highest order differ- 
ential operator, the linear analogue of (1.1) is singular as we described above. 



1 This second condition may be added with no loss of generality and slightly simplifies the 
analysis. See for instance [37]. 
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With this in mind we expected that the full nonlinear model might also be sin- 
gular. As developed below, v may be found exactly from a single observation 
over a finite time window which suggests that this singular structure is preserved 
in this nonlinear case. 

The starting point of our analysis, the derivation of an estimator for v, 
follows methods already developed for the linear case (see references mentioned 
above). We project (1.1) down to a finite dimensional space, and for each TV we 
arrive at a system of the form 

dU N + (uAU N + P N B{U))dt = P N adW, U(0) = U , 

where Pn is the projection operator on the finite dimensional space generated 
by the first TV Fourier eigenvalues of the Stokes operator. We then formally 
compute the MLEs associated to these systems, and take them as an ansatz for 
our estimators. In the course of the analysis we introduce an additional degree 
of freedom, a parameter a, which we may carefully tune to compensate for the 
nonlinear term. We arrive finally at the following three classes of estimators: 

f*(A 1+2a U N , dU N } + fi(A 1+2a U N , P N B{U))dt 

J T \A 1 + a U N \ 2 dt 
J^{A 1+2a U N , dU N ) + J^{A 1+2a U N , P N B{U N ))dt 

fi\Ai+"U«\*dt (1.5) 

N 

E^ +2a K(T)-u 2 (0)-T\- 2 r) 

k=l 

N _ 

fc=i 

Here Uk ■= {U, $&) represents the fcth (generalized) Fourier mode of the solution 
U. 

The main result in this work establishes the following properties for the 
proposed estimators: 

Theorem 1.1. Suppose that U = U(uj) is a single sample path solution of (1.1 ), 
(1.2) or (1.1), (1-3) observed on a finite interval of time [0,T]. Assume that 
(1.1) is forced with a white noise process of the form (1-4) where 7 > 1, and 
suppose that a > 7 — 1. Then, given a suitably regular initial vector field Uq, 

(i) the functionals z/jv, vn, >^n defined by (1.5) are weakly consistent estima- 
tors of the parameter v, i.e. 

lim = lim vn — hm = v 

N->oo N->oo N—>-oo 

in probability. 



VN = - 



v N = - 



v N = - 



2 In the case (1.3) we assume, for technical reasons, an upper bound on 7, 7 < 1 + 1/4 as 
well. See below. 
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(ii) if we assume further that a > 7 — 1/2, then vn is asymptotically normal 
with rate N i.e. 

N(v N -v) — ► 77 

(converges in distribution) where n is a mean zero, normally distributed 
random variable. 

While we are able to prove the strongest convergence results for vm, this 
estimator is intractable numerically and even analytically. This is because vn 
depends on all of the Fourier modes of the solution in a highly nonlinear fashion. 
At the other extreme is £jv which is much more straightforward to compute but 
is expected to have a slower rate of convergence to the actual parameter v. The 
estimator i>n is a compromise between the two extremes since it depends only 
on the knowledge of the first N eigenmodes but retains some of the complex 
structure of the nonlinear term. Although at the present time we are not able 
to prove this, we expect that vn has a faster rate of convergence than v^- We 
conjecture, in Section 4.3, that is also asymptotically normal with the same 
variance and rate of convergence as v^. Given the explicit formulas for the 
estimators, (1.5), all these questions, including the effect of the free parameter 
a on the rate of convergence, can be studied by means of numerical simulations, 
which the authors plan to undertake in a separate forthcoming paper. 

While the form of the proposed estimators and the general statements of 
the main results in this work are similar to previous works in the linear case, 
fundamental new difficulties arise which require one to take a novel approach 
for the analysis. This is of course due to the complex structure of the nonlinear 
term appearing in (1.1) which couples, in an intricate way, all of the modes 
Uk = (U, In contrast to the linear case, we lose for example any explicit 
spectral information about the elements Uk- This coupling also means that the 
Uk are not expected to be independent. 

To overcome these difficulties the analysis relics on a careful decomposition 
of the solution U = U + R. Here U satisfies a linear system where the modes 
are independent. Crucially, a complete spectral picture is obtainable for U . On 
the other hand, R, while depending in a complicated way on the full solution 
U is more regular in comparison to U. This is because R is not directly forced 
by the noise terms adW. For this point the analysis, particularly in the case of 
bounded domains, requires a delicate treatment of the nonlinear term. 

Due to these technical issues, we were able to establish asymptotic normality 
only for vn- It is interesting that vn is a consistent estimator for v and it is the 
same as the MLE of the corresponding linear equation (the stochastic Stokes 
equation). This effect can be explained as follows: since the nonlinear term 
B(U) is in some sense 'lower order' it fails to destroy the information about v; 
v remains observable in a similar manner to the linear case. 

The exposition of the paper is organized as follows: In Section 2 we lay 
the theoretical foundations for this work reviewing the relevant mathematical 
theory for the stochastic Navier-Stokes equations. We establish some crucial 
spectral information concerning the linear system associated to (1.1). We also 
recall in this section some particular variants on the law of large numbers and 
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the central limit theorem. Section 3 sketches the derivation of the estimators 
v, v, vn. We conclude the section with a strict formulation of the main results. 
The proof of the main theorem is carried out in Section 4 in a series of modular 
substeps. We first study the regularity of the 'residual' R that appears after we 
'subtract off' the noise term appearing in (1.1) via the linear Stokes equation. 
As an immediate application we arc able to determine some precise rates for 
the denominators appearing in the estimators (1.5) . Using these rates we 
successively analyze the consistency of the estimators. The final subsection 
treats the question of asymptotic normality with the help of a central limit 
theorem for martingales. 

2 Mathematical Setting of the Problem 

We begin by recalling the mathematical background for the stochastic Navier- 
Stokes Equations and then review some general results from probability theory 
that will be used in the sequel. 

2.1 The Stochastic Navier-Stokes Equation 

We first describe how (1.1) is recast as an infinite dimensional stochastic evolu- 
tion equation of the form 

dU + (uAU + B(U))dt = adW , 

/ \ (2.1) 

U(0) = u . y 1 

The basic functions spaces are designed to capture both the boundary conditions 
and the divergence free nature of the flow. 

We first consider the spaces associated with a Dirichlet boundary condition 
(1.3). Let H := {U G L 2 {V) 2 : V ■ U = 0, U ■ n = 0}, where n is the outer 
pointing unit normal to dV. H is endowed as a Hilbert space with the L 2 
inner product (U b , 0) = J v U b U i dx and associated norm \U\ = (U, U) 1/2 . The 
Leray-Hopf projector, Pjj, is defined as the orthogonal projection of L 2 (T>) d 
onto H. We next take V := {U G Hq(V) 2 : V • U = 0} and endow this space 
with the inner product ((C/ b ,C/ J )) = f M ^U b ■ V0dM. Due to the Dirichlet 
boundary condition, (1.3), the Poincare inequality \U\ < c\\U\\ holds for U G V 
justifying this definition. 

The definitions for H and V are slightly different for the case of periodic 
boundary conditions (1.2). We take T> = [— L/2, L/2} 2 and define the spaces 
L 2 er (T)) 2 , Hp er (T>) 2 to be the families of vector fields U = U(x) which are L 
periodic in each direction and which belong respectively to L 2 (0) 2 and i/ 1 (0) 2 
for every open bounded set 0CK 2 . We now define 

H = (u eL 2 per {V) 2 :V -U = 0, J t/(x)dx = oj , 

and 

V =\u e H^ er {V) 2 :V-U = 0,J U(x)dx = J . 
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H and V are endowed with the norms | • | and || • || as above. Note that we impose 
the mean zero condition for H and V so that the Poincare inequality holds. As 
mentioned in the introduction, there is no loss of generality in imposing this 
extra assumption. See, e.g. [37]. 

The linear portion of (1.1) is captured in the Stokes operator A = — Ph A, 
which is an unbounded operator from H to H with the domain D(A) = H 2 (A4)n 
V. Since A is self adjoint, with a compact inverse A" 1 : H — > D(A), we may 
apply the standard theory of compact, symmetric operators to guarantee the 
existence of an orthonormal basis {&k}k>i for H of cigenfunctions of A with 
the associated eigenvalues {Xk}k>o forming an unbounded, increasing, sequence. 
Moreover, 

Afe^Aifc, (2.2) 

where the notation a„ « b n means that limn-,.^ a n /b n = 1. Also, we will write 
in ~ b n when there exists a finite, nonzero constant c such that lim„_ ) . 00 a n /b n — 
c. For more details about asymptotical behavior of {Xk}k>i see for instance 
[1, 29] for the no-slip case (1.3), and [6] for the spatially periodic case (1.2). 
Define Hn = Span{$i, . . . , &n}, and take Pn to be the projection from H onto 
this space. We let Qn '■= I — Pn- 

The analysis below relies extensively on the fractional powers of A. Given 
a > 0, take D(A a ) = {U g H : Y.k X T\ u k? < °°}, wh erc u k = (U,$ k )- On 
this set we may define A a according to A a U = J2 k X^Uk^k, for U = J2 k Uk&k- 
Classically we have the generalized Poincare and inverse Poincare estimates 

\A°»P N U\ < X%- ai \A^P N U\, \A°*Q N U\ < -J^\A a *Q N U\, (2.3) 

X N 

for any oi\ < ct2- 

We next describe the stochastic terms in (1.1). Fix a stochastic basis S := 
(0, J 7 , {J-t}t>0)P, {Wk}k>i), that is a filtered probability space with {Wk}k>i 
a sequence of independent standard Brownian motions relative to filtration 
{J~t}t>o- In order to avoid unnecessary complications below we may assume 
that Tt is complete and right continuous (see [8] for more details). Writing for- 
mally W = J2k>o ®kWk, W may be viewed as a cylindrical Brownian motion 
on H. 

We briefly recall the classical formalism for infinite-dimensional Wiener pro- 
cess as in [8], [34]. Consider the collection of Hilbcrt-Schmidt operators mapping 
H into D(A /3 ), /3 > 0. We denote this family by L 2 {H 1 D{A fs )). Throughout 
this work we assume that a, understood as an operator, has the form 

a$ fe = Ap$ fc . (2.4) 

We will write adW{t) = J2k>i X k 7 $kdW k (t), t > 0. One may check that, 
for every e > 0, a g L 2 (H, D^" 1 / 2 "')). In particular, given the standing 
assumption that 7 > 1, we have a g L-2(H, D{A 1 / 2 )). 



Parameter Estimation for the Stochastic Navier-Stokes Equations 



9 



2.2 The Stochastic Stokes Equation and Limit Theorems 

We next consider the linear system associated to (2.1), which we write in the 
abstract form: 

dU + vAUdt = ^A~ 7 $ fc cW fc , U(Q) = U . (2.5) 

k 

For the purposes here this system can be analyzed as 2D stochastic heat equation 
driven by an additive cylindrical Brownian motion (for general results we refer 
readers to [8, 35].) 

Let us denote by u k , k > 1, the Fourier coefficients of the solution U with 
respect to the system {&k}k in H, i.e. u k = {U, $fe), k > 1. By (2.5), we 
note that each Fourier mode u k represents a one dimensional stable Ornstein- 
Uhlcnbcck process with dynamics 

du k + v\ k u k dt = \~~<dW k , u k (0) = ii ak , k > 1. (2.6) 

It follows from (2.6) that 

u k (t) = u k (0)e- vXkt + A" 7 / e-^-^dWkit), fc > 1, i > 0. (2.7) 



In what follows we will use the following auxiliary results about asymptotics 
of the first moments of the Fourier modes u k , k > 1 (see also Theorem 2.1 in 
[26].) 

Lemma 2.1. Suppose that U is a solution of (2.5) and let U N := P N U, N > 1. 
(i) Assume that 7' < 7 and that E| A 7 '~ 1/2 £/ | 2 < 00. Then 

U £ L 2 (n; LL([0, 00); D(A 7 '))) n L 2 (n- C([0, 00); D(A-''-^ 2 ))) . (2.8) 
(ii) Suppose that Uq = 0, then: 

r T Ti -(l+2 7 ) ^.-(1+27) 

E / u 2 k dt » ±±h a (2 .9) 

./n 2^ 2^ 1 ; 



Var 



T 







A; (3+47) ~ k -^\ (2.10) 



fra) Moreover, for /3 > 7, 
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Proof. The first item is classical and may, for example, be justified with a 
Galerkin Scheme or other suitable techniques from the general theory of ex- 
istence and uniqueness of the solutions for stochastic parabolic equations. Sec 
e.g. [8, 35]. Using (2.7), (ii) follows by direct computations of the corresponding 
moments, and for the final item we deduce 

El \A f3 U N \ 2 dt = E f \J2\ p k u k <S> k \ 2 dt = J2\fE[ u\dt 
Jo Jo fc=l fc=l J o 

_ J_ ,2,9-1-27 ^ 1 A l iV 

~ 2v^ k ~ 2 7 2/3 - 2-y ' 



where we have made use of (ii), (2.2) in conjunction with 

N N 1+a 

fe=i 

The proof is complete. □ 

We hnally recall some particular versions of the Law of Large Numbers (LLN) 
and the Central Limit Theorem (CLT) which are used to prove consistency and 
asymptotic normality of the class of estimators given by (1.5). 

Lemma 2.2 (The Law of Large Numbers). Let £„, n > 1, be a sequence of 
random variables and b n , n > 1, an increasing sequence of positive numbers 
such that lim n ^oo b n = +oo, and 



^ Varg„ 

Z^-^2— <0 °- ( 2 - 13 ) 



n=l 

(i) If we assume that the random variables £„, n > 1, are independent then 



E (6 - Mk 

k=l 
n— >oo fo. 



lim ^— i — = a.s. 



(ii) If we suppose only that £„, n > 1, are merely uncorrelated random vari- 
ables, then 

n 

lim — — = 0, (2.14) 

n->oo b n 

in probability. 

Proof. See, for example, Shiryaev [36, Theorem IV. 3. 2] for the proof of (i). Part 
two, similar to the proof of Weak LLN, follows from the Markov inequality. For 
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a fixed e > 0, and for all pairs m < n, we have 



fe=i 



> e 



\ 



\ 1 

^ E I - ) ^ ^ E Va < 



6 2 



k=l 



m 1 n 

^^E Var ^ + ^E 



Var^ 



k=l 



k—rn 



1 m 1 

^^E Var ^ + ^E 



fe=i 



*2 



Since 6„ — > oo, (2.14) follows. 



□ 



The following central limit theorem is a special case of a more general result 
for martingales; see, for instance [24, Theorem 5.5.4(11)]. 

Lemma 2.3 (CLT for Stochastic Integrals). LetS = (fi, J",P, {J"t}t>o, {W fc } fc >i) 
6e a stochastic basis. Suppose that o~k £ L 2 (f2; L 2 ([0, T])) is a sequence of real 
valued predictable processes such that 



N 



lim 



E Jo °idt 
fc=i 



1 m Probability. 



fc=i 



Ther 



N 



E Jo 



A' 



fc=i 



1/2 



converges in distribution to a standard normal random variable as N — > oo. 



2.3 The Nonlinear Term 

The nonlinear term appearing in (2.1) is given by B(U, [/") := Ph{(U ■ V)C/*) = 
PHdfj^UjdjUt), which is defined for £7 6 V" and £7* e D(A). Note that, 
for brevity of notation, we will often write B(U) for B(U, U) as for example in 
(2.1). We have the following properties of B: 

Lemma 2.4. 

(i) B is bilinear and continuous from V x V into V and from V x D(A) into 
H . For U, £7" £ V, B satisfies the cancelation property 



(B{U,U t ),U t )=0. 



(2.15) 
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If U, £/", U b are elements in V, then 

{B(U,U s ),U b ) < c|[/| 1/2 ||C/|| 1/2 ||f/ B |||[/ b | 1/2 ||^H 1/2 - (2.16) 
On the other hand if U G V , £ D(A), and U b G H, then we have: 

(B(U,U*),U b ) 



< c 



|C/|i/2||£/||i/2||C/ll||i/a|^ll|i/a|^b| 
|[/|i/2|^|i/2 ii^lliii^bi, 



(2.17) 



(ii) In the case of either periodic, (1.2) or Dirichlet, (1.3) boundary conditions, 
B(U) £ D(A f} ) for every < f3 < 1/4 and every U £ D(A). Moreover, 
for such values of (3, 



\A?B{U)\ 2 < c\\U\\ 2 \AUf 



(2.18) 



(Hi) In the case of periodic boundary conditions (1.2), B(U, U*) £ D(A°) when- 
ever $ > 1/2, U £ D(AP), U" G D{A fi+1 / 2 ), and for such U, U%, 



\A»B(U,U*)\ 2 < c\A^U\ 2 \A^ 2 U\ 2 . 



(2.19) 



Proof. The properties outlined in (i) and (iii) are classical; see, for instance, 
[38], or [6, Lemma 10.4] for (2.19). 

The properties in (ii) are established via interpolation and the equivalence 
of certain fractional order spaces, see [14]. Since [14] emphasized the case of 
spatial dimension 3, for the sake completeness, we briefly recall the arguments. 

For any element U £ D(A), standard estimates imply that 

\B(U)\ 2 < c\\Uf\AU\, 

\\B{U)f H1{M)2 <c\\U\\\AU\ 3 . 

Let V = H n i? 1 (I?) 2 and, for s £ (0,1) we define the interpolation spaces 
V s = [V,H]i_ 8 . Sec [23] for the general theory. In [14], it is established that 
D(A@) = V2/3, (3 < 1/4 in Dirichlet case (1.3). 3 Note that V does not incorporate 
boundary conditions and so B(U) £ V, for U £ D(A). In consequence, for any 
such U £ D(A) and allowed values of0</3<l/4we have, by interpolation 

\Af>B(U)\* = \B(U)\l w < (|B(t0r- a/, ||5(C0ll^) 3 
<c(||J7|| 3 |A[/|) 1 - 2/3 (||[/|||AJ7| 3 ) 2 ' 3 
< c(\\U\\ 3 \AU\) 1/2 (\\U\\\AU\ 3 ) 1/2 < c||J7|| 2 |A[/| 2 . 

Combining these observations gives (ii), completing the proof. □ 

3 In the periodic case, (1.2), V = V so that D(A@) = V 2 p, as a direct consequence of the 
fact that D(A 1 / 2 ) = V. 
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Remark 2.5. When we consider the case (1.3) it is not true in general that 
B(U) € D(A^), even for U G D(A^ +1 ' 2 ), [3 > 1/4. This is due to the fact that 
while the Leray projector Pg is continuous on H m (T>), m > 1, we do not expect 
that P H maps H™(V) into H™(V). See [38] and also [14]. For this reason we 
may not expect an inequality like (2.19) for such Dirichlet boundary conditions. 
As such, (2.18) relies on a delicate analysis of small fractional order space where 
the boundary is not present; see [14, 23]. 

2.4 Existence, Uniqueness and Higher Regularity 

With these mathematical formalities in place we now define precisely (2.1), in 
the usual time integrated sense and recall some now well established existence, 
uniqueness and regularity results for these equations. Note that for this work 
the solutions we consider correspond to so called 'strong solutions' in the deter- 
ministic setting (see [38]). In the context of stochastic analysis, since we may 
suppose that the stochastic basis S is fixed in advance, we may say that the solu- 
tions considered are 'strong' (or less confusingly 'pathwisc') in the probabilistic 
sense as well. 

Theorem 2.6. 

(i) Suppose that we impose either (1.2) or (1.3) and assume that Uq G V, 
a G L 2 (H,V). Then there exists a unique, H -valued, J- t -adapted process 



with the equality understood in H . 
(ii) In the case of periodic boundary conditions (1.2) if ft > 1/2 so that a G 



L 2 {H,D(AP)), Uo G D(A?), then 

U G LL«0, oo), D(A^)) n C([0, oo), D{A?)). (2.21) 



Remark 2.7. 

(i) As noted above, when a is defined via (2.4), a G L 2 {H, V) whenever 7 > 1. 
Indeed we have a G L 2 (H, D{A^)) for every /3 < 7 — 1/2. 

(ii) We suspect that higher regularity similar to Theorem 2.6, (ii) may be 
established in the case of Dirichlet boundary conditions, (1.3). However 
since (2.19) docs not apply (sec Remark 2.5) a different proof than outlined 
here is needed. 



U with 



U G LL([0, 00); D{A)) n C([0, 00); V) a.s. 



(2.20) 



and so that for each t > 0, 
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Proof. The well-posedness of (2.1) has been studied by many authors as dis- 
cussed in the introduction. Since we are considering the case of an additive noise 
the proof is close to the deterministic case after we perform a suitable change of 
variables. For completeness, we briefly recall some of the formal arguments and 
note that the computations may be rigorously justified with a suitable Galcrkin 
scheme. Consider hrst the linear system (2.5) with initial condition U(0) = Uq. 
As in Lemma 2.1 above, we have that U in Lf oc ([0, oo); D(A)) D C([0, oo); V) 
(or in Lf oc ([0, oo), D(A^ +1 / 2 )) n C([0, oo), D{AP)), under the conditions of item 
(ii)). We now consider the shifted variable U = U — U , which satisfies 

^-+vAU + B(U + U) = U(0)=0. (2.22) 
dt 

The estimates that lead to (2.21) are standard. We hrst multiply (2.22) by 
U, integrate over the domain and use (2.15), (2.16), (2.8) to infer that U E 
Lf oc ([0, oo); V) (~1 Lf£ c ([0, oo); H). With this regularity in hand we next multiply 
(2.22) by AU and apply (2.17), (2.8) in order to conclude (2.20). 
For ft > 1/2 we multiply (2.22) by A 2l3 Uand infer 

d| ^f |2 + 2v\A fi+1 ' 2 U\ 2 - 2{A' 3 B(U + U), A?U) = 0. (2.23) 
Since (3 > 1/2, we may apply (2.19) and estimate 

d ^L + 2„\A^ 2 U\ 2 
dt 

<c\A l3 (U + U)\\A fs+1 ^{U + U)\\A p U\ 
<c(\A^U\ 2 + \A p tf\ 2 )\A p U\ 2 + u\A 0+l,2 U\ 2 + v\A 0+r/2 U\ 2 . 



Rearranging, 

d\A^U\ 2 
dt 



■ 2v\AP +1 / 2 U\ 2 < c(\Af>U\ 2 + \A?U\ 2 )\APU\ 2 + v \A^/ 2 Uf 



Observe that, due to the Gronwall Lemma, if U £ Lf oc ([Q, oo); D(A^)) then we 
infer that U £ L? oc ([0, oo); L>(A /3+1 / 2 ))nL^ c ([0, oo); D{A P )). The desired result 
therefore follows from an inductive argument on /3 starting with the base case 
assumption j3 € [1/2, 1) which is satisfied as a consequence of (2.20). □ 

3 Estimators for v\ Heuristic Derivation and the 
Main Results 



In this section we sketch the heuristic derivations of the estimators based on a 
particular version of the Girsanov Theorem. We then restate, now in precise 
terms, the main results of this paper. 
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As before we denote by U N the projection of the solution U of the original 
equation (2.1) onto Hn = PnH = M. . Note that U satisfies the following 
finite dimensional system: 

dU N = -{vAU N +^ N )dt + P N cxdW, U N {0) = U^, (3.1) 

where tp N ' (t) := Pn(B(U)). To obtain an initial guess of the form of the estima- 
tor for the parameter v, we treat ip N as an external known quantity, independent 
of v and view (3.1) as a stochastic equation evolving in M. N . Let us denote by 
P^' T the probability measure in C([0, T]; R w ) generated by U . Formally, we 
compute the Radon-Nikodym derivative or Likelihood Ratio dP^ ,T /dF^ T (see 
e.g. [25, Section 7.6.4]) 

dP?> T (U N ) 



dl 



N,T ~ eX P 



(- j n (v-v Q )(AU N )'G 2 dU N (t) 

(v 2 -v 2 )(AU N )'G 2 AU N dt 
{v -v n ){AU N )'G 2 ^ N dt) 1 



where G := (Pno')^ 1 = diagfuf 1 , . . . , er^ 1 ] = diag[A7, • • • , Ajy] and v' denotes 
the transpose of the vector v G Mr. By maximizing the Likelihood Ratio with 
respect to the parameter of interest v, we may compute the (formal) Maximum 
Likelihood Estimator (MLE) of the parameter v. A direct computation 
yields 

J T {AU N )'G 2 dU N + £(AU N YG 2 P N (B(U))dt 

vn = t- ■ (o.^ ) 

J^(AU N )'G 2 AU N dt 

As expected, vn is a valid estimator and in fact one can show that it is a 
consistent estimator of the true parameter v. This consistency makes essential 
use of the fact that the denominator [AU N )' G 2 AU N dt diverges to infinity 
as N t oo. Sec Lemma 4.3 below. With this in mind, wc introduce a slight 
modification to the MLE (3.2), and propose the following class of estimators 

J T (A^U N ,dU N )+J Q T (A^U N ,P N B(U))dt 

vn = -fp , (0.6) 

J c , | A 1 + a U N \ 2 dt 

where a is a free parameter with a range specified later on. Note that this 
formulation appears in the functional language developed above and is derived 
using that the action of G 2 on Hm is equivalent to A 21 ' . Also we observe that 
vn is a particular case of vn with a = 7. 

While the estimator vn has desirable theoretical properties, it also assumes 
that Pn(B(U)) is computable, which could be quiet a difficult task. Since our 
goal is to provide estimators that can be eventually implemented in practice 
(evaluated numerically), wc propose two further classes of estimators. One class 
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is naturally derived from (3.3) by approximating Pn{B(U)) with Pn(B(U )) 

/ r (i4 1 + 2a l7 JV ) dU N ) + £(A^U N , P N B(U N ))dt 

vn = t . (3.4) 

\A 1 + a U N \ 2 dt 

Note that now depends only on the first N Fourier modes. However, even in 
this case the expression for PnB(U n ) is very complicated due to the nontrivial 
coupling of the modes. See e.g. [10]. It turns out, as shown rigorously below 
(see Proposition 4.6), that the second term appearing in (3.3), 

fi(AWU N ,P N B(U))dt 

kn ■= -— ~r , (3.5) 

£\A l +<*U N \Ht 

is of lower order and tends to zero, as N — > oo. Hence we get the following 
consistent estimators of the parameter v 

N 1 2a T 

Jv(A 1+2a U N ,dU N ) £ Afe UkdUk 

£>n = 



fe=i 



(3.6) 



N 

Z\l +2a (ul(T)-ul(0)-T\ft) 
fc=l 

k=l 

Clearly this last estimator is easiest to compute numerically. On the other hand 
it may lack the speed of convergence of the first two. 

We conclude this section with the main result of this paper: 

Theorem 3.1. Suppose that U solves (2.1) with either (1.2) or (1.3) in the 

sense of and under the conditions imposed by Theorem 2.6. Assume that 7 > 1 
and in the case (1-3), additionally that 7 < 1 + 1/4. Also, assume that Uq € 
D{AP), for some ft > 7 - 1/2. 

(i) If a > 7 — 1, then v^, vn and i>n as given by (3.3), (3.4), and (3.6) are 
weakly consistent estimators of the parameter v, i.e. 

lim = lim vjy = lim On = v 

N-yoo JV-Voo TV— s-oo 

in probability. 

(ii) If a > 7 — 1/2, then is asymptotically normal with rate N, i.e. 

N(v N -v)^n, (3.7) 
where rj is Gaussian random variable with mean zero and variance 

2t/(ct-7+l) 2 
AiT(q- 7 +1/2) ' 
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4 Proof of the Main Theorem 

We establish the proof of the Theorem 1.1 in a series of propositions. As men- 
tioned in the introduction, we do not have precise spectral information about 
Fourier coefficients Uk = (U, k > 1, in contrast to linear case (see Section 
2.2). To overcome this, we proceed by decomposing the solution into a linear 
and a nonlinear part, U = U + R. We assume that U is the solution of the linear 
stochastic Stokes equation (2.5) with U(0) = 0. The residual R must therefore 
satisfy, 

d t R + vAR = —B(U), R(0) = R o . (4.1) 

First, we study the regularity properties of R and show that R is slightly 
smoother than U. Subsequently we make crucial use of this extra regularity 
and establish the consistency of the proposed estimators by showing that second 
term in (3.3) converges to zero. The final section treats the asymptotic normality 
using CLT introduced in Section 2.2. 

Remark 4.1. For simplicity and clarity of presentation we shall assume a more 
regular initial condition Uo G D(A 1 ) in contrast to the statement of Theo- 
rem 3.1. The more general case when we assume merely that U$ € D(A") for 
some /? > 7 — 1/2 may be treated by writing U = U + R + S, where U satisfies 
(2.5) with Uq = 0, R satisfies (4.1), this time with i?o = 0, and finally 5 is the 
solution of d t S + vAS = 0, with 5(0) = U . 



4.1 Regularity Properties for the Residual 

Proposition 4.2. Suppose that (3 > 0, 7 > 1 and that R solves (4.1) with U 
the solution of (2.1) corresponding to an initial condition Uo S D{A 1 ^ 2+ ^). 

(i) If U and R satisfy Dirichlet boundary conditions (1.3), and (3 < 1/4, then 
for every T > we have 

sup [AW+PRl 2 + ( T \A 1+ ?R\ 2 < 00. (4.2) 
te[o,T] Jo 

Moreover, for an increasing sequence of stopping times T n with t„ j" 00, 

E ( sup lA 1 / 2 ^! 2 + [ T " \A 1+p R\ 2 ) < 00. (4.3) 
\te[o,T„] Jo J 

(ii) In the case that both U and R satisfy periodic boundary conditions (1.2) 
and we assume (3 < 7 — 1/2 the same conclusions hold. 

Proof. As above in Theorem 2.6 the computations given here may be rigoursly 
justified via Galerkin approximations. Multiplying (4.1) by A 1+2 °R, integrating 
and using the symmetry of the powers of A we infer 

1 ^1^+1/2^2 + ^+1^2 = -( A P B (U), A 0+1 R). (4.4) 
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For the case of a bounded domain, (1.3), we infer from (2.18) and Theorem 2.6, 
(2.20) that, 

/ \A p B(U)\ 2 dt<c \\Uf\AU\ 2 dt < oo a.s. 
Jo Jo 

By integrating (4.4) in time and making standard estimates with Young's In- 
equality, (4.2) now follows in this case. 

In the case of the periodic domain we estimate \A B(U)\ differently. Dehne 
/3' = max{/3, 7/2} so that 1/2 < /3' < 7 — 1/2. By applying the higher regularity 
estimates (2.19) we find that 

(A p B(U), A^ +1 R) I <\A p B(U) | \A P+1 R\ 
<c\A p 'B(U)\\A fi+1 R\ 
<c\A^U\\A^' +1 / 2 U\\A f3+1 R\ 

<c\a^'u\ 2 \a^ +1 ' 2 u\ 2 + ^\AP +1 R\ 2 . 

Due to Theorem 2.6, (ii), we have, for any T > 0, that 

[ T \A^U\ 2 \A^' +l ' 2 U\ 2 dt< 00 a.s., 
Jo 

and (4.2) follows once again. 

For the stopping times t„, we define 

r n := inf jsup \\U\\ 2 + J \AU\ 2 dt' > n\ 

when (1.3) is assumed and 

r„ := inf jsup \A 15 ' U\ 2 + \A> 3 ' +1 ' 2 U\ 2 dt' > n\ 

for (1.2). In either case it is clear that {T n } n >i is increasing. Moreover, in the 
case (1.3), since P(r„ < T) = P(su Pt ,< T ||J7|| 2 + J T \AU\ 2 dt' > n), it follows 
from (2.20) and the fact that r n is increasing, that lim n _> r n = 00 a.s. Arguing 
in the same manner for the case (1.2), the proof is complete. □ 

Remark. Comparing Proposition 4.2, (4.2) with Lemma 2.1, (2.8), (2.11) we 
see that R has been shown to be just shy of a derivative more regular than U. 
More precisely we have that, for any e > 0, 

U G L 2 (n; L 2 OC ([0, 00); D(A^))), U # L 2 (n ; L 2 OC ([0, 00); D{A^))), 

while on the other hand, 

i?ei 2 (O;LL([0,»);^(A 7+1/2 " E ))). 
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As an immediate application of these properties of the residual R we have 
the following result: 

Lemma 4.3. Suppose that U and U are the solutions of (2.1) and (2.5) respec- 
tive^ For both (1.2), and (1.3) we suppose that 7 > 1, U(0) = U G D(A 7 ) 4 
and [To = 0. Additionally, in the case (1.3), we assume that 7 < 1 + 1/4. Then, 
for any a > 7 — 1 , 

Urn Aj! 1 = 1 (4.5) 

^E/ T \A 1 + a U N \ 2 dt 

with probability one. 
Proof. Note that 

\A 1+a U N \ 2 <\A 1+a U N \ 2 + lA^JZ"! 2 + 2\A 1+a U N \\A 1+a R N \, 

| A l+ a[/ JV|2 >| A l+a^JV|2 + \ A l+a R N\2 _ 2 \A 1+a U N \\A 1+a R N \ , 

and therefore (4.5) follows once we have shown that 

Ej^\A 1 + a U N \ 2 

and that 



Jo ■ — ^ ~ 



2 ,_ ^ rT 



E/ lA^l/^l 



2 



a.s. (4.6) 



For the first item, I± , we apply the law of large number (LLN), Lemma 2.2, 

T n 

with £„ := A 2a+2 / u 2 n (t)dt and 6„ := £ E [6J- Notice that, due to, (2.2) and 

fe=i 

(2-9) 

n n 

fo « ~ E x T +2x k 1-27 ~ E fc 2Q " 27+1 . (4.7) 



fe=i fc=i 



Given the assumptions a > 7 — 1, we have that lim„_ ! . 00 b n = 00. Moreover, 
combining (4.7) with (2.2), (2.10), (2.12) 

V X^ii „ V A «" +4 ~ 3 ~ 47 
52 ~ 



fc=i 



A 4a-4 7 +l j 



~ X! /\2q-2 7 +2n 2 ~ n 3 < °°' 

n>l W > n>l 

Thus, by the LLN we conclude that lini/v_ i . 00 if = 1 with probability one. 



4 At the cost of further evaluations, this condition may be weakened to the conditions 
imposed in Theorem 3.1. This applies both here and below for Propositions 4.4, 4.6. See 
Remark 4.1. 
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Since 1 + a > 7, by (2.11) we infer 

E [ T \A 1+a U N \dt~N 2a -^+ 2 . 
Jo 



Pick any a' G (7— 1, min{a, 1/4}), in the case (1.3), or any a' G (7— 1, minja, 7— 
1/2}) under the assumption (1.2). By applying (4.2) for R established in Propo- 
sition 4.2, we have in both cases that 

T 

\A 1+a ' R\ 2 dt < 00 a.s. 

Combining these observations and making use of (2.3), we have 

J Q T \A^R N \ 2 dt \ 2 ^- a '^^ A i +a > R N ?dt fi\ A i+ a ' R?dt 

2 — C pj2a-2f+2 — C J\f2a-2f+2 — ]\[2a'-2y+2 

Due to the restrictions on the choice of a', we have that 2a' — 2-f + 2 > 0, and 
hence Jjj — > 0, as N — > oo, with probability one. The proof is complete. □ 

4.2 Consistency of the Estimators 

Using the dynamics of U N , i.e. substituting (3.1) into (3.3), we get the following 
representation for the estimator vn'- 



C{A 1+2a U N ,P N adW) 

vn = v — y — = 

£\A 1 + a U N \ 2 dt 



j T (A^+ 2 ^u N ,j:% =1 <P k dw k ) 

= v — 

Similarly, we deduce 



(4.8) 



\A 1 + a U N \ 2 dt 



(4.9) 



S,{A^ 2 ^u\Y. N k =i^dw k ) 

VN = V 7p 

Jo \A 1 + a U N \ 2 dt 
; tf(A^ 2a U N ,P N B(U)-P N B(U N ))dt 
J T \A 1 + a U N \ 2 dt 

Note that vn = un — kn, with kn defined by (3.5). Thus, 

Jo (A 1+2a ~ J U N , Y^k=i ®kdW k ) 

v N = v - n N 1 = . (4.10) 

J^ \A 1 + a U N \ 2 dt 

With the above representations for the estimators, the consistency will follow 
if we show that each stochastic term on the right-hand side of (4.8), (4.9), (4.10) 
converges to zero. 
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Proposition 4.4. Assume the conditions and notations from Lemma 4-3. Then, 
(i) for every 8\ < min{2 + 2a — 27, 1}, 

. C {A 1+2a ^U N <$> k dW k ) , 
w->oo \A 1 + a U N \ 2 dt 

(ii) whenever 62 < min{2 + 2a — 2j, 3/2} m the case (1.2), or whenever 62 < 
min{2 + 2a — 27, 5/4 + 1 — 7} in the case (1.3), we have 



K^T^?*"*-* (4.12) 



jS 2 Jo 



in probability. 



Proof. Due to Lemma 4.3, (4.5) and (2.11) the desired result follows once we 
show that each of sequences 



J 

and 



1 . = 5»(A 1+2a ^U N ,Y. N k= ^dW k ) Ht^' 1 Si n k dW k 

N ■ ,2+2a-27-5i < 2+2a-27-<Si 

N A N 



T 2 ._ Io(A 1+2a ^R N ,E k=1 ^dW k ) 

N ■ . 2+2a-27-<5 2 

A N 

converge to zero as N — > 00. 

For the first term, J*., define £ k := X 1 + 2a ^ f Q T Uk dW k and b„ := 
Under the given conditions, lim„_>.oo b„ = 00. With the Ito Isometry and (2.9), 
we have 

Varfe] = E[|£] ^ A 2+4a-2 7A -(l+2 7 ) = A l+4a-4^ 



Thus, 



hi ~ x 4+4q-47-2(Si Z^ \3-2i5i ~ Z_^ ^3-2.51 00 " 

n>l An n>l An n>l 



Note that under the given conditions 5\ < 1. This justifies the assertion that the 
final sum is finite. We conclude, by the LLN, Lemma 2.2 that liniAr^oo = 0. 

We turn to Jjy. Let r k := (R,§ k ), k > 1, and for any stopping time r we 
define 

C k :=X k +2a ^ [ T r k dW k . 



Note that the random variables (1, k > 1, are uncorrelated. Similarly to the 
above arguments, we let b„ := ^2+20-27-53 an( j 0Dserve that this sequence is 
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increasing and unbounded. Up to any stopping time r such that Var££ < oo, 
we have 



E 52 fe — E \4+4a-4 7 -2<5 2 E /_ T ^ t 
k>l k k>l A k 



2+4a-27 



Va^ 2+252 E f r\dt (4.13) 
Jo 



fe>i 

; / \A~<- l+&2 R\ 2 dt 



Note that under initial assumptions, in the case of a bounded domain, (1.3), 
7— 1 + 82 < 5/4, and in the periodic case, (1.2), we have 7 — l + <5 2 < 7 + 1/2. In 
either case, by taking r„ as in Proposition 4.2, we infer from (4.13) with (4.3) 
that, for every n, 



E 

/c>i 



Var C r r " 



6 : 



< 00. 



By applying Lemma 2.2, we conclude that, for each n fixed, 

hm — ttt^ — n — 1 = 111 Probability. 

JV^-oo ,2+2q-2 7 -,5 2 J 

N 

Since r n is increasing, f2 = U n {r„ > T} is a set of full measure, and a simple 
estimate yields that 

l im — . 2+2a - 2 T-t - = ° inProbabilit y- 

iV 



The proof is complete. □ 

Corollary 4.5. Putting the admissible values Si = 82 = in (4.11), (4.12), and 
taking into account that U N = U N + R N , we conclude 



lim ^ = '^ fc=1 =i =0 in Probability. 



/„ [il 1 + a I7 w | 2 dt 



Thus, by representation (4.8) we have that vn is weakly consistent estimator of 
the true parameter v. 

We turn next to the 'nonlinear terms' appearing in (4.9). 

Proposition 4.6. Assume the conditions and notations imposed for Lemma 4-3 
above. We suppose that 8 € [0, min{5/4 — 7, a — 7+ 1}) in the case (1.3) or that 
8 S [0, min{l/2, a — 7 + 1}) when we assume ( 1.2). Then 

, ^{A 1+2a U N ,P N B{U))dt 
lim N S Jo N „ V ;/ = a.s. 4.14 

N^oo jT \A 1 + a U N \ 2 dt 
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Proof. By the Cauchy- Schwartz inequality, 



fi(A l +*°U N ,P N B{U))dt 
^\A l + a U N \ 2 dt 



fi\A*P N B{U)\ 2 dt 
J T \A 1 +°>U N \Ht 



1/2 



< 



(4.15) 



Due to Lemma 4.3, (4.5) and (2.11) it is therefore sufficient to show that 



We begin with the boundary conditions (1.2) and consider two possibilities 
corresponding to different values of a. First suppose that a < 7 — 1/2, so that 
5- (a- 7 + 1) < 0. Pick any /3 € (max{a, 1/2}, 7- 1/2). Making use of (2.19) 
and then applying Theorem 2.6, (ii), we observe that 



and (4.16) follows. 

Now suppose that a > 7 — 1/2. In this case we pick an element a' E 
(max{<5 + 7 — 1, 1/2}, 7 — 1/2). Note that, by assumption, 5 < 1/2 so that this 
interval is non-trivial. Clearly a' < a and we apply (2.3) and again (2.19) in 
order to estimate 



As above we find that the quanity on the right hand side is finite due to The- 
orem 2.6, (ii). Noting that S — (a — 7 + 1) + a — a' < 0, we infer that (4.16) 
holds true. 

The case of Dirichlet boundary conditions (1.3) is addressed in a similar 
manner. When a < 1/4 we directly apply (2.18) to infer (4.16). When a > 1/4 
we pick any a' G (5 + 7—1,1/4). Noting that the conditions on S ensure 
that this interval is nontrivial and that a' < a, we apply (2.3) and (2.18), in 
a similar manner to (4.17) and infer (4.16) for this case too. The proof is now 
complete. □ 

Corollary 4.7. In similar manner one can establish the same results as above 
for PnB(U n ). In particular, for 5 = we have 




(4.16) 





lim 

N^oo 



lim 

N—HX 



^(A^ a U N ,P N B(U))dt 
Jo \A 1+a U N \ 2 dt 
^{A 1+2a U N ,P N B{U N ))dt 
J T \A 1 + a U N \ 2 dt 



a.s. 



a.s. 



(4.18) 



(4.19) 
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Taking into account the above equalities and the representations (4.9) and (4.10) 
we have that and vn are consistent estimators of v. 

4.3 Asymptotic Normality 

We finally address the asymptotic normality of vjq and prove second part of 
Theorem 3.1. Using the representation (4.8) for vn. Lemma 4.3, and (4.5) we 
see that is suffices to establish that 

where ri is normal random variable with mean zero and variance, . 2 ^"~ 7 ^V/on 

' ' Ai I (a — 7+1/2) 

and that 

lim 7V^- Tf, _ — =0 in Probability. (4.21) 

N ^°° E/ T \A 1 +»U N \ 2 dt 

We establish (4.20) with the aid of Lemma 2.3. Let a k := A^ +2a ~ 7 Ufc, and 
£fc : = Jo a l dt > k>l. Notice that, due to (2.9), 

E[&] ~ A 2+4q " 27 A; (1+27) = \\+^ , 

Varfe] ~ A 4+8«-47 A -(3+4 7 ) = A l+8a-8 7 _ 

Define b n := Xwb=i ^4k- Under the given assumptions, 1 + 4a — 4~f < —1, so 
that by (2.12) we have that b n ~ A 2+4q_47 . We infer that b n is increasing and 
unbounded. Moreover, 

k=l °* k=1 

and therefore by LLN, Lemma 2.2, we conclude 

JV 

lim — -F = 1 a.s. 

JV-s-oo N 

fc=i 

Consequently, by Lemma 2.3 with a k defined above, we have 

■b. f'^^-^^UArtO,!). (4.22, 



Noting that both 1 + a > 7, 1 + 2a — 7 > 7 we may apply (2.11) and infer 
(EjT\A^^U»\ 2 dt) 1/2 r^r a _ 1 + l 1 



E/olA^^Pdi V AiT ^a- 7 + 1/2 iV 



(4.23) 



Parameter Estimation for the Stochastic Navier-Stokes Equations 



25 



By combining (4.22) and (4.23) we obtain (4.20). 

The other condition, (4.21), follows directly from Lemma 4.4, (4.12) with 
the admissable value of §2 = 1. This completes the proof of Theorem 3.1, (ii). 

Remark 4.8. Notice that the accuracy of the estimators z/jv, as measured by 
the variance from asymptotic normality, depends on a, 7 and T. As one would 
expect, as T gets larger and more information is revealed, the quality of the 
estimator improves (i.e. the variance decreases). This suggests that one may 
show that each of the classes of estimators considered above are consistent in 
the large time asymptotics regime, i.e. when we fix N and send T — > 00. As 
with the spectral method we have developed in this work establishing the long 
time asymptotics is more complicated in comparison to the linear case and will 
be addressed in future work. 

Also, we note that (a — 7 + 1)/ \J a. — 7 + 1/2, as a function of a on the 
domain a > 7 — 1/2 reaches its minimum at a = 7. Thus, for fixed T and 7, 
the smallest asymptotic variance for the estimator vn corresponds to a = 7. 
Observe that when a = 7 the estimator vn reduces to the formal MLE (3.3), 
which in some sense is the optimal estimator in this class of estimators. 

Remark 4.9. We want to emphasize that we still believe, that asymptotic 
normality properties similar to (3.7) also hold true for the estimators 0^. How- 
ever, for a rigorous proof one needs to show, for example, that Proposition 4.6 
holds true for some 5 > 1, and with P N B{U) replaced by P N B{U) - P N B{U N ). 
Intuitively it is clear that the difference PnB(U) — PnB(U n ) will make the 
convergence to zero in (4.14) faster, allowing for a larger 8 comparative to those 
from the terms PpfB(U) and PnB(U n ) considered individually. We further 
believe that the quality of these estimators vjq may be optimized in terms of 
the free parameter a. Although, up to the present time, we remain unable to 
establish such quantitative results about the asymptotic normality of the esti- 
mators vjq, we plan to study these questions at least by means of numerical 
simulations in forthcoming work. 
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